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Abstract. The non-autonomous chiral model equation for an m x m matrix function 
on a two-dimensional space appears in particular in general relativity, where for m = 2 
a certain reduction of it determines stationary, axially symmetric solutions of Einstein's 
vacuum equations, and for m = 3 solutions of the Einstein-Maxwell equations. Using 
a very simple and general result of the bidifferential calculus approach to integrable partial 
differential and difference equations, we generate a large class of exact solutions of this 
chiral model. The solutions are parametrized by a set of matrices, the size of which can 
be arbitrarily large. The matrices are subject to a Sylvester equation that has to be solved 
and gcnerically admits a unique solution. By imposing the aforementioned reductions on 
the matrix data, we recover the Ernst potentials of multi-Kcrr-NUT and multi-Dcmiariski- 
Newman metrics. 
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1 Introduction 



The bidifferential calculus framework allows to elaborate solution generating methods for a wide 
class of nonlinear "integrable" partial differential or difference equations (PDDEs) to a conside- 
rable extent on a universal level, i.e. resolved from specific examples. It takes advantage of 
the simple rules underlying the calculus of differential forms (on a manifold), but allows for 
a generalization of the latter, which is partly inspired by noncommutative geometry. For a brief 
account of the basic structures and some results we refer to pQ (also see the references therein) , 
but all what is needed for the present work is provided in Section [2j In this framework we 
explore the non-autonomous chiral model equation 

(pfeO,+(w 1 ) P = ° ( L1 ) 

for an mxm matrix g, where p > and z are independent real variables, and a subscript indicates 
a corresponding partial derivative. It apparently first appeared, supplemented by certain reduc- 
tion conditions (see Section [5]) , as the central part of the stationary axially symmetric Einstein 
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vacuum (m = 2) and Einstein-Maxwell (m = 3) equations (see in particular [21 El SI El El d EJ ) - 
For m > 3 this equation is met in higher-dimensional gravity, with a correspondingly enlarged 
number of Killing vector fields (see e.g. [21 HOI HU H21 HH1 [III [HI HH HZl [13 HH])- A version 
of the above equation also arises as the cylindrically symmetric case of the (2 + l)-dimensional 
principal chiral model [20] and as a special case of the stationary Landau-Lifshitz equation for 
an isotropic two-dimensional ferromagnet [21j . 

The first construction of "multi-soliton" solutions of (jl.ip has been carried out by Belinski 
and Zakharov [51 E] (also see [7]) using the "dressing method" 0. Here JHJ| is expressed as 
the integrability condition of a linear system, which depends on a (spectral) parameter and 
involves derivatives with respect to the latter. Another approach is based on a linear system 
that depends on a variable spectral parameter, i.e. a parameter that depends on the variables p 
and z [3]. In Appendix [B] we show that both linear systems arise from a universal linear system 
(see Section [2]) in the bidifferential calculus framework (also see [31] for a relation between the 
two linear systems). 

In the present work, we concentrate on a surprisingly simple general solution generating result 
in the bidifferential calculus framework, which has already been successfully applied in various 
other cases of integrable (soliton) equations [H [32| [33"1 [33] to generate multi-soliton families. 
In order to make it applicable to the non-autonomous chiral model, a slight generalization is 
required, however (see Section [3] and Appendix [A]) . Section U] then elaborates it for the m x m 
non-autonomous chiral model. We obtain solutions parametrized by four matrices. Two of them 
arise as solutions of an n x n matrix version of the quadratic equation for pole trajectories that 
first appeared in the solution generating method of Belinski and Zakharov [5[ El [7J. It then 
remains to solve a Sylvester equation, where two more matrices enter, which are constant of 
size m x n, respectively n x m. Since n can be arbitrarily large, we obtain an infinite family of 
solutions. The Sylvester equation is easily solved if the first two matrices are chosen diagonal, 
and in this case one recovers "multi-soliton" solutions. Additional solutions are obtained if the 
two n x n matrices are non-diagonal. In this case it is more difficult to solve the Sylvester 
equation, though a not very restrictive spectrum condition ensures the existence of a unique 
solution. Except for an example in Section EJ we will not elaborate this case further in this 
work. 

Section [5] addresses reductions, in particular to the Ernst equation of general relativity. It 
turns out that the "multi-soliton" solutions of the stationary, axially symmetric Einstein vac- 
uum and Einstein-Maxwell equations are indeed in the generated class of solutions of the non- 
autonomous chiral model. We thus obtain a new representation of these solutions. It has the 
property that the superposition of two (or more) "solitons" (e.g. black holes) simply corre- 
sponds to block-diagonal composition of the matrix data parametrizing the constituents. This 
puts a new perspective on an old result about one of the most important integrable equations 
in physics. 

We would like to stress that the solutions of the non-autonomous chiral model and the 
Ernst equation(s), (re)derived in this work, originate from a universal result that also generates 
multi-soliton solutions of various other integrable equations in a non-iterative way. The crucial 
step is to find a "bidifferential calculus formulation" of the respective equation. This may be 
regarded as a generalization of the problem of formulating the equation as a reduction of the 
selfdual Yang-Mills equation. Indeed, in the case under consideration, it is of great help that 
an embedding of the non-autonomous chiral model in the (m x m) selfdual Yang-Mills equation 
is known [351 EH [37J [381 EH HO] , arid a bidifferential calculus formulation is then obtained from 

Further constructions of multi-soliton solutions, in the context of general relativity, were presented in par- 
ticular by Alekseev |22l [23] , Neugebauer [24] , Kramer and Neugebauer |25l 126] , Korotkin [27] (limiting cases of 
finite-gap solutions), Manko et al. |28II29] , and Masuda et al. [30] (using Hirota's method). Also see the references 
in [8] Section 34] for other solution generating techniques. 
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that of the selfdual Yang-Mills equation pQ, see Section [U Once this is at hand, the remaining 
computations are rather straightforward. 

Section contains some concluding remarks. 

2 Preliminaries 

Basic definitions. A graded algebra is an associative algebra f2 over C with a direct sum 
decomposition Q = © r>0 ^ r into a subalgebra A := 0° and ^4-bimodules fT, such that O r ft s C 
Q r "M. bidiff eren tial calculus (or bidifferential graded algebra) is a unital graded algebra Q 
equipped with two (C-linear) graded derivations d, d : $7 — > Q, of degree one (hence dU r C Q r+1 , 
dQ r C f2 r+1 ), with the properties 

= Vk£C, where d K := d - «d, (2.1) 

and the graded Leibniz rule 

d K (xx) = (d K xW + (-l) r xd K x', 

for all x £ ^ r an d x' G ^- This means that d and d both satisfy the graded Leibniz rule. In 
Section [3] we consider a more narrow class of graded algebras. A bidifferential calculus within 
this class is then specified in Section UJ 

Dressing a bidifferential calculus. Let (fi, d, d) be a bidifferential calculus. Replacing d K 
in ([21]) by 

D K := d — A — «d, 

with a 1-form A (i.e. an element of O 1 ), the resulting condition T)\ = (for all k£C) can be 
expressed as 

dA = and dA - AA = 0. (2.2) 

If these equations are equivalent to a PDDE or a system of PDDEs for a set of functions, we 
say we have a bidifferential calculus formulation for it. This requires that A depends on these 
functions and the derivations d, d involve differential or difference operators. There are several 
ways to reduce the two equations (I2.2p to a single one. Here we only consider two of them. 
1. We can solve the first of (I2.2p by setting 



A = d(f). 

This converts the second of (|2.2[) into 

dd^ = d<M</>. (2-3) 

This equation is obviously invariant under i— > acf)a~ l + /? with an invertible a G A satisfying 
da = da = 0, and j3 £ A satisfying d/3 = 0. 

2. Alternatively, the second of equations ()2.2p can be solved by setting 

A = (d<7)<r\ 

and the first equation then reads 

d((d 5 ) 5 - 1 ) = 0. (2.4) 
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This equation has the (independent left and right handed, i.e. chiral) symmetry 

g ^ ag/3, (2.5) 

where a G A is d-constanlH and (3 G A is d-constant, and both have to be invertible. Since 

d[(d 5 - 1 ) 5 ] = -dig-'dg) = g- 1 [(ddg)g- 1 - (dg^g^g = g^d^dg^g, (2.6) 

g solves fj2.4|) iff g" 1 solves (|2.4p with d and d exchanged. In our central example, the non- 
autonomous chiral model, g i— > g" 1 becomes a symmetry. 

Linear system. The compatibility condition of the linear equation 

dX = (dX)P + AX (2.7) 



is 



= d 2 X = (dX)[(dP)P - dP] + (dA - A 2 )X - (dA)XP. 
If P satisfies 

dP = (dP)P, (2.8) 
this reduces to 

(dA - A 2 )X = (dA)XP. (2.9) 



For the above choices of A, this implies the respective PDDE. Hence (|2.7p is the source of 

a corresponding Lax pair, also see Appendix [Bl 

As a consequence of (j2.8|) . P solves (j2.3|) and, if P is invertible, also (|2.4p , 

Miura transformation. If a pair (</>, <?) solves the Miura transformation equation 

(d^ 1 = d</> (2.10) 



(cf. [I]), it follows (as an integrability condition) that (j) solves (|2.3p and g solves (|2.4p . We 
note that (|2.10p is just the linear equation (|2.7p if we identify A = d<f>, X = g and set P = 0. 
If we have chosen a bidifferential calculus and a reduction condition such that (12.3P becomes 
equivalent to some PDDE, this does not necessarily mean that also (j2.4j) is equivalent to some 
"ordinary" PDDE. But for the central example of this work, the non-autonomous chiral model, 
such a mismatch does not occur. In fact, in Section [3] we will actually present a solution 
generating method for (I2.10p . 



3 A solution generating method 

Let /\(C N ) denote the exterior (Grassmann) algebra of the vector space and Mat(m,n, B) 
the set ofmxn matrices with entries in some unital algebra B. We choose A as the algebra of 
all finite-dimensional matrices (with entries in B), where the product of two matrices is defined 
to be zero if the sizes of the two matrices do not match, and assume that f2 = A <8> /\(C ) is 
supplied with the structure of a bidifferential calculus. In the following, I = I m and I = I n 
denote the m x m, respectively n x n, identity matrix. 

2 Although not evident, we need not require da = in addition. 
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Proposition 3.1. Let P,R,X £ Mat(n, n,B) be invertible solutions of 
dP = {dP)P, dR = RdR, 

dX = (dX)P - (dR)X, XP - RX = VU, (3.1) 
with d- and d-constant U E Mat(m, n,B), V £ M&t(n,m,B). Then 

cf) = UX- 1 V, g = I + U(RX)- 1 V (3.2) 

solve the Miura transformation equation ()2.10|) . and thus (j2.3j) . respectively ()2.4|) . 

Proof. Using the last three of (|3.ip we obtain 

d(RX)' 1 = -X^fiXX- 1 + R^dRlR- 1 = -X~ 1 (dX)X~ 1 (XP)(RX)- 1 
= (dX- 1 )[J + VU(RX)- 1 }. 

Multiplication by U from the left and by V from the right, and using dl = 0, leads to 

dg = UidX-^Vg = (dfig. 

Hence (f> and g solve the Miura transformation equation (|2.10p . We did not use the first of (|3.ip , 
but it arises as an integrability condition: = d 2 X = (dX)[(dP)P - dP]. U 

Remark 3.1. The third of (|3.1|) . which has the form of the linear equation (|2.7|) . is almost 
a consequence of the fourth, which is a Sylvester equation. Indeed, as a consequence of the 
Sylvester equation we have 

= d{RX - XP + VU) = {dR)X + RdX - {dX)P - XdP 
= R[dX + (dR)X] - [dX + XdP]P 

= R[dX + (dR)X - (dX)P] - [dX + (dR)X - (dX)P]P + d(RX - XP)P, 

where the last term vanishes. If P and R are sufficiently independent, this implies that the 
third of (|3.ip is satisfied. In particular, this holds if B is the algebra of complex functions of 
some variables and if P and R have no eigenvalue in common. 

Appendix [X] explains how Proposition 13.11 arises from a theorem that has been applied in 
previous work to generate soliton solutions of several integrable PDDEs. 

4 The non-autonomous chiral model 

The PDE defining the non-autonomous chiral model can be obtained as a reduction of the 
self-dual Yang-Mills (sdYM) equation (see e.g. [351 ESI EH [381 EHl SO] ) . In an analogous way, 
a bidifferential calculus for the non-autonomous chiral model can be derived from a bidifferential 
calculus for the sdYM equation (also see [H]). In coordinates p, z, 9, where p > 0, it is given by 

df = -hd + e e (f P - /T7*)C 2 , df = e~\f p + p~ l fe)Qi + f&. (4-1) 

Here e.g. f z denotes the partial derivative of a function / (of the three coordinates) with respect 
to z, and d, £2 is a basis of /\ 1 (C 2 ). d and d extend to matrices of functions and moreover to 
U = A®/\(C 2 ) with A = Mat (m, m, C), treating (j, £2 as constants. The coordinate 9 is needed 
to have the properties of a bidifferential calculus, but we are finally interested in equations for 
objects that do not depend on it. 
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A (matrix- valued) function is d-constant (d-constant) iff it is ^-independent and only depends 
on the variables 9, p through the combination pe s (respectively pe~ e ). It is d- and d-constant 
iff it is constant, i.e. independent of z, 9, p. 

For an m x m matrix-valued function g, (|2.4p takes the form 

{p9 z g~ 1 ) z + {pg P g~ l ) p - {g P g~ l ) e + {9eg^) p - p~ l {g e g~ l ) d = o. 

Restricting g by setting 
g = e ce g 

with any constant c and ^-independent g, for the latter we obtain the non- autonomous chiral 
model equation^] 

{p~gz~g- l ) z + {p~g P ~g- l ) p = ®- (4.2) 

In Section 14.11 we derive a family of exact solutions by application of Proposition 13.11 In 
Appendix [B] we recover two familiar linear systems (Lax pairs) for this equation. 
Miura transformation. Evaluating (|2.3p with 



where 4> is ^-independent, we obtain 

4>zz + 4> PP + p~ x 4> P = [4> P + p' 1 ^, 4> z ] , (4.3) 

which is related to the non- autonomous chiral model by the Miura transformation 

4>z = -g P g~ l - cp^i, g z g~ l = 4> p + p" 1 ^. 

Symmetries. (|4.2h is invariant under each of the following transformations, and thus, more 
generally, any combination of them. 

(1) g i— > ag(3, with any invertible constant m x m matrices a and (3 (cf. (|2.5p ). 

(2) g I—)- p c g with any constant c. 

(3) ~g ^ g' 1 (also see flZS) ). 

(4) g i — y g^ , where ' indicates Hermitian conjugation. 

We note that g i— > (^^)~ 1 is a fairly obvious symmetry. With its help, (4) follows immediately 
from (3). 

4.1 A family of exact solutions 

Let us first consider the equation dP = (dP)P, which is the first of (|3,ip . Using the above 
bidifferential calculus, it takes the form 

P Z P = -e~ 9 {P p + p- l P e ) , P z = e 9 (P p - p^Pg) P. 

Writing 

P = e~ e P, 

and assuming that P does not depend on 9, this translates to 

P P ~ P~ l P = -PzP, Pz = {Pp + P~ 1 P)P. (4.4) 

The proof of the following result is provided in Appendix O 

3 Changing the sign of the first term in the expression for d/ in (|4.1|l . we obtain a minus sign between the 
two terms on the left hand side of (|4.2|) . This hyperbolic version of the chiral model shows up, in particular, in 
the reduction of the Einstein vacuum equations with two spacelike commuting Killing vector fields, describing 
gravitational plane waves [7]- Our further analysis can be adapted to this case. 
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Lemma 4.1. The following holds. 



(1) 



If P and I + P are invertible, the system (|4.4p implies 



P 2 - 2p~ 1 {zl + B)P -1 = 0, 



(4.5) 



with a constant matrix B. 

(2) Let I + P 2 be invertible and P p , P z commute with P. If P satisfies (|4.5[) . then P 
solves (|4.4|) . 

Remark 4.1. If P is diagonal, then (|4.5[) becomes the set of quadratic equations (2.11) in [6] 
(or (1.67) in |7]), which determine the "pole trajectories" in the framework of Belinski and 
Zakharov. In our approach, there are more solutions since P need not be diagonal. 

Remark 4.2. The conditions [P p ,P] = [P Z ,P] = in part (2) of the lemma are satisfied in 
particular if the spectrum spec(B) is simple, i.e. if the eigenvalues of B are all distinct, since 
then the solutions of (|4.5p are functions of p, z and the matrix B (and thus P p and P z commute 
with P) [42]. But this would be unnecessarily restrictive, see Section [4.21 

Remark 4.3. Under the assumption that I + P 2 is invertible, (|4.4|) implies [P p , P] = [P Z , P]=0, 
also see (jC.l[) . For the bidifferential calculus under consideration, dP = (dP)P is therefore 
equivalent to dP = PdP. The latter is one of our equations for R in Proposition 13. 11 Setting 

R = e~ e R, 

with R ^-independent, invertible P and R both have to solve (|4.5p . 

The third of (|3.ip becomes 

X p + p-^Xe = -X Z P + R Z X, X z = (X p - p- l X 9 )P - (R p + p^R)X. 

Assuming that U and V are ^-independent, and recalling the 0-dependence of <f>, the formula 
for <fi in (|3.2p requires X = e e X with ^-independent X. Hence 

X p + p~ l X = -X Z P + R Z X, X z = (X p -p' 1 X)P-(R p + p~ 1 R)X. (4.6) 
The last of (]3.ip becomes the ^-independent Sylvester equation 

XP - RX = VU. (4.7) 
Now Proposition 13.11 implies the following. 

Proposition 4.1. Let n x n matrices P and R be solutions of (j4.5|) (with a matrix B, respec- 
tively B 1 ), with the properties that they commute with their derivatives w.r.t. p and z, and that 

~ 2 ~ 2 ~ 

I + P and I + R are invertible. Furthermore, let spec(P) n spec(i?) = and X an invertible 
solution of the Sylvester equation (|4.7p with constant m x n, respectively n x m, matrices U 
and V . Then 



with any constant invertible m x m matrizQ go, solves the non- autonomous chiral model equa- 
tion (EHj) . 




(4.8) 



4 Here go represents the freedom of chiral transformations. 
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Proof. As a consequence of the spectrum condition, a solution X of the Sylvester equation (|4.7p 
exists and is unique. The further assumptions for P and R are those of Lemma 14.11 part (2). 
Furthermore, (|4.6p is a consequence of (|4.7p if the spectrum condition holds (also see Remark l3.ip . 
Now our assertion follows from Proposition 13.11 and the preceding calculations. ■ 



Remark 4.4. The determinant of (|4.8p is obtained via Sylvester's theorem, 

det(s) = det (/ + UiRX^V) det(g ) = det (/ + VUiRX)" 1 ) det(g ) 

= det(RX + VU) det(-RX)- 1 det (50) = det(XP) det(flX)" 1 det (go) 
det(P) , . , 

where we used the Sylvester equation (|4.7p and assumed that it has an invertible solution. 

Remark 4.5. As an obvious consequence of (|4,7p . U and V enter g given by (|4.8p only modulo 
an arbitrary scalar factor different from zero. We also note that a transformation 

P^T~ 1 PT 1 , R^T~ 1 RT 2 , Ut->UTi, V ^ T~ l V , X^T^XT^ 

with constant invertible nx n matrices Ti, T2, leaves (|4.5p . (I4.6p . (|4.7j) and (I4.8P invariant. As 
a consequence, without restriction of generality, we can assume that the matrix B in (|4.5p . and 
the corresponding matrix related to R, both have Jordan normal form. 

Example 4.1. Let P and R be diagonal, i.e. 

P = (Pi5ij) , R=(ri5ij). 

If they have no eigenvalue in common, then (|4.7p has a unique solution given by the Cauchy-like 
matrix 



Xij 



(VU)ij 



Pj - n 

It remains to solve (|4.5p (choosing B diagonal), which yields 



Pi = p- 1 (z + b i +j i y/(z + b i ) 2 + f?), r i = p - 1 [z + b' i +j' i J{z + b' i y + f?), (4.9) 



with constants 6j, b[ and G {±1}- Since we assume that {p{\ n {r^} = 0, the assumptions 
of Proposition 14.11 are satisfied. It follows that, with the above data, (|4.8j) solves the non- 
autonomous chiral model equation. 

The case where P or R is non-diagonal is exploited in the next subsection. But Example 14. II 
will be sufficient to understand most of Section [5j 

4.2 More about the family of solutions 

Introducing matrices A and L via 

A = (zI + B) 2 + p 2 I, P = p~ l (L + zI + B), 
(14. 5p translates into 
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According to Remark 14,51 we can take B in Jordan normal form, 

B = block-diag(.B ni , . . . , B n J. 
Let us first consider the case where B is a single r x r Jordan block, 

/0 1 ••• o\ 

1 



B r = bl r + N r , N r 



Then we have 

A = x 2 (I r + M r ), 

where 

M r = x- 2 [2(z + b)N r + N*], 
and thus 

r-l 



1 



L = x(I r + M r ) 1 / 2 = xJ2( 1/ k 2 )M 



by use of the generalized binomial expansion formula, noting that M r r = as a consequence of 
N r r = 0. Hence we obtain the following solution of (j4.5|) . 



Pr = P~ l (zlr + B r + r ( X f ) M^j 



(4.10) 



which is an upper triangular Toeplitz matrix. In particular, we have 
P l = p- 1 [z + b + t], 
P 2 = p- x [z + b + x\ 

y u -L / 

\p 2 x~ 3 (z + b + r)" p 



p 1 [z + b + x] 













r- 1 


e 


1 







/i 


r- 1 





1 














1 

^ 2 r _3 (2; + 6 + r) 

1 




-l l 



^{z + b) p 2 x- 5 {z + b + x)~ 1 \ 
lp 2 x~ 3 (z + b + x)~ 1 



These matrices are obviously nested and, from one to the next, only the entry in the right upper 
corner is new. 

For the above Jordan normal form of B, solutions of (|4.5p are now given b}0 
P = block-diag(P ni , . . . , P n J, 



5 For P with simple spectrum, every solution of (|4.5[) has this form. Otherwise, i.e. when there are two Jordan 
blocks with the same eigenvalue, (|4.5|l has additional solutions, see [32]. But they are further constrained by (|4.4p . 
We will not consider such solutions in this work. 
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where the blocks typically involve different constants replacing b, i.e. different eigenvalues of B. 
Since P p and P z obviously commute with P, and since I+P is invertibk^l for p > 0, Lemma l4.ll 
part (2), ensures that P solves (|4.4p . If P has the above form, and R a similar form, and if P 
and R have disjoint spectra, it remains to solve the Sylvester equatiorjjl (|4.7|) in order that fj4.8j) 
yields solutions of the non-autonomous chiral model equation. This leads to a plethora of exact 
solutions. We postpone an example to Section[5l where additional conditions considerably reduce 
the freedom we have here, see Example 15.21 

5 Reductions of the non- autonomous chiral model 
to Ernst equations 

According to Section 01 a particular involutive symmetry of the non- autonomous chiral mo- 
del (|4.2p is given by g \— > 7 (<?^) _1 7, where 7 is a constant matrix with 

7 t = 7) 7 2 = J. 

(|4.2p therefore admits the generalized unitarity reduction g^jg = 7, which means that g belongs 
to the unitary group £7(m; 7 )H Another reduction, associated with an involutive symmetry, is 
= g. Imposing both reductions simultaneously, amounts to setting 

~g ] = 9, d~g) 2 = i- (5.1) 

Writing 

~g = l(I-2V), 
translates these conditions into 

7 pt 7 = V, V 2 = V. 

In particular, V is a projector. If we require in addition that rank("P) = 1, which for a projector 
is equivalent to tr('P) = 1 |431 Fact 5.8.1], the following parametrization of g can be achieved 
(also see e.g. m EE2 EE3 E]), 

~g = J-2—, (5.2) 

where v is an m-component vector with v^v 7^ 0. This parametrization is invariant under 
v 1 — y cv with a nowhere vanishing function c, so that the first component of v can be set to 1 in 
the generic case where it is different from zero. If 7 has signature m—1, (|5.2p is a parametrization 
of the symmetric space SU(m-l, l)/S(U(m - 1) x U(l)) jSJ 051 SB] . The condition tr(V) = 1 
corresponds to 

tr( 7 g) = m — 2. (5-3) 

We also note that det(<?) = — det( 7 ). The following result, which we prove in Appendix[Cl shows 
how the reduction conditions (|5.ip and (|5.3p can be implemented on the family of solutions of 
the non-autonomous chiral model obtained via Proposition 14.11 



~ 2 



r-l 

3 Note that I r + PI = (1 + p~ 2 (z + b + r) 2 )J + Yl QkN 1 ; with functions g k . 

k-l 

7 Under the stated conditions the Sylvester equation possesses a unique solution and a vast literature exists to 
express it. 

8 If 7 has p positive and q negative eigenvalues, this is commonly denoted U(p,q). 
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Proposition 5.1. Let X solve the Sylvester equation (|4.7p . where P, R, U, V, satisj 

(TP) 2 = -I, (TR) 2 = -I, g o7 U = UT, TV = Vg 0l , (5.4) 

with an n x n matrix T and a constant m x m matrix go satisfying 

r 2 = I, (sot) 2 = (5-5) 

Furthermore, let spec(P) n spec(-R) = 0. 

(1) g given by (|4.8|) satisfies 

(rrg) 2 = I and tr( 7 5)=tr( 75o )-2tr(r). (5.6) 

(2) If moreover the relations 

R^ = TPT, U^ = Vg , gl=g , T^ = T (5.7) 

hold, then g given by (|4. 8[) is Hermitian. 

Remark 5.1. Let the matrix data (Pi, Ri, Ui, Vi, Tj) satisfy T 2 = I ni and 

(r 4 P 4 ) 2 = (T t R t ) 2 = -I ni , go jU i = U i T i , T i V i = V i g j. 

Set P = block-diag(Pi,...,Pjv), R = block-diag(fli , . . . , R N ), V = block-diag(ri, . . . , F N ), 
and 



U = (U 1 ,...,U N ), V 



\v N/ 



Then we have T 2 = I and (I5.4p holds. If spec(P) n spec(i?) = 0, the corresponding Sylvester 
equation has a unique solution X. According to part (1) of Proposition 15. 1| g given by (|4.8p 
satisfies the reduction conditions ()5.6|) . This is a way to superpose solutions from the class 
obtained in Sectional preserving the constraints (|5.4|) . We simply block-diagonally compose the 
matrix data associated with the constituents. In an obvious way, this method can be extended 
to part (2) of Proposition 15.11 

Remark 5.2. Let n = 2N and 

where P and -R are invertible block-diagonal N x N matrices, composed of blocks of the 
form (|4.10p . Then we have (TP) 2 = —I n and (TR) 2 = —I. Choosing 7 and go such that 
(gol) 2 = I) the conditions gojU = UT and TV = Vgoj are solved by 

U=(U ig 01 U), V=(_.^ o7 ), 
where U and V are arbitrary constant m x N, respectively N x m matrices. Writing 

x = (x tc x zp\ 

" \ Z RXP J ' 



These conditions are motivated by the structure of P, R, U, V found in Example 15.1 
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reduces the 2N x 2N Sylvester equation (|4.7p to the two N x N Sylvester equations 

XP - RX = VU, ZP + R^Z = -iVgo-fU. (5.8) 



If X and Z are invertibel, then RX is invertiblq 10 !. Proposition 15. 11 part (1), implies that 
with the above matrix data satisfies ('yg) 2 = I and ti(yg) = tr(7po)- With a suitable choice 
of 7 and go we can achieve that (|5.3p holds. To fulfil the remaining Hermiticity condition, one 
possibility is via part (2) of Proposition 15.11 See also Examples 15.21 and 15.71 Such solutions can 
be superposed in the way described in Remark 15. 11 

In the special case where R = rljsr, the solutions of the Sylvester equations (|5.8p are X = 
VU(P - rlTv)" 1 and Z = -iVg jU(P + r -1 /^) -1 . These expressions are not invertible 
if N > m, so in this particular case our solution formula only works for N < m (also see 
Example E2}. 



Remark 5.3. The following observation in particular underlies the Harrison transformation 
2] which we consider in Example 15.51 below. Let H be an m x m-matrix that satisfies 



H^H = 7. (5.9) 
If g satisfies 

iig) 2 = i, a ] = 9, 

then also 

if = H~gH\ 

and we have tr(7^') = tr(7^). If g has the form (14. 8p with Hermitian g$ and (jgo) 2 = I, and 
if H also satisfies 

Hg H^ = g , (5.10) 
then the effect of the transformation g i— >• g' amounts to the replacement 

U h-> U' = HU, V ^V' = VH~ 1 , 
which leaves the Sylvester equation (|4.7p invariant. 

5.1 Solutions of the Ernst equation of general relativity 

We choose m = 2 and 



7 = 

and write 

v = 



i 
-i 



with a complex function £ and its complex conjugate E. Then (|5.2p takes the form 



9 



i Us-s) 



8 + 8\^{S-£) ££ 

10 (RX)^ 1 can be computed as a 2 x 2 block matrix. The problem of evaluating the original expression for g 
that involves 2N x 2N matrices then reduces to that of evaluating only N x N matrix expressions. 
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so that 

s= 1 - i ~ 921 . 

(|4.2p now becomes the Ernst equation 

(Re£){d 2 p + p~ l d p + d 2 z )8 = (8 p f + (£ z ) 2 , 

where e.g. d p denotes the partial derivative with respect to p. This equation determines solutions 
of the stationary axially symmetric Einstein vacuum equations. The following statements are 
easily verified. 



1. Excluding g = ±7, the reduction conditions (15. ip are equivalent to 

g real, det(g) = 1, 512 = 521- (5.11) 



()5.3[) is then automatically satisfied. 

2. For real g, the second of the reduction conditions (|5.ip implies the first. As a consequence, 
Proposition 15.11 part (1), already generates solutions of the Ernst equation. 

We will use these observations in the following examples. 

Example 5.1 (Kerr-NUT). For the solution of the non-autonomous chiral model given in 
Example 14.11 with n = 2, we have (also see Remark 14. 4 p 

det(g) = ^^det(go), 
r±r 2 



with pi, Ti given by (|4.9p . Choosing 
9o = h, 

so that det(po) = 1> the second of the reduction conditions (|5.1ip is solved by setting 

1 1 

P2 = , r 2 = , 

Pi n 

noting that —1/pi is given by the expression for pi with jj exchanged by — jj. We shall write p, r 
instead of p±, r\. With U = (uij), V = (%•), the remaining constraint = g 2 i is solved bjo 

u 2 2 = -U1XU12/U21, v 2 2 = -viiv 2 i/v 12 - 
In the following we assume that un and vu are different from zero and write 

U21 = UUll, V\2 = vvu. 

Then uu, u±2, «n, V21 drop out of g. Without restriction of generality, we can therefore choose 
them as u\\ = 1, U12 = —u, v\\ = 1 and V21 = —v, hence 

Then U and V commute with 7. g is reaZ in particular if either of the following conditions is 
fulfilled. 



Another solution is U22 = Ui2W2i/«ii,«22 = Wi2«2i/«ii. But this leads to p = go- 



14 



A. Dimakis, N. Kanning and F. Miiller-Hoissen 



(1) p, r real (which means b, b' real) and u, v real. 

(2) f = — - (which means b' = b and j = —f G {±1}) and v 

The Ernst potential takes the form 



(1 + uv) 



p+r 



1 pr- 
\Pr- 



(u — i)(v 



By a shift of the origin of the coordinate z, we can arrange in both cases that 

p = p- 1 ( z + b + jt+) , r = p- 1 (z - b + j't_) , 
t± := V(z±6) 2 +p 2 , j,/e{±l}, 

where 6 € R in case (1) and 6 £ iR in case (2). Using 



p + r 



1 



p — r 2b 
and introducing 

, 1 + uv 



(jx + + fx. 



pr — 1 



1 



pr + 1 26 



(jt+ - ft-), 



(5.12) 



(5.13) 



u — v 



I = jb- 



uv 



m 



we obtain 



i|(t++ii / r_)-2(m + iQ 



i§(t4 +j/t_) + 2(m + i0' 



(5.14) 



Setting jf = — 1, the cases (1) and (2) simply distinguish the non-extreme and the hyperextreme 
Kerr-NUT space-times (see e.g. |5Uj). The constants satisfy m 2 + I 2 — o 2 = b 2 . 

Example 5.2. Let n = 2N. With the choices made in Remark 15.21 Proposition 15.11 part (1), 
implies that g, given by (|4.8p with go = I 2 , satisfies {^g) 2 = I 2 and tr(7<?) = 0. Choosing 
all parameters real, it follows that g = I 2 + U(RX)~ 1 V determines a solution of the Ernst 
equation, provided that X is invertibld* 2 !. For N = 1, we are back to the preceding example. 
For N = 2 let, for example, 



where p 



(p px 1 





° \ 


p 











-p - 


1 (pt)- 1 


\0 






/ 


i 0^ 







i 




-i 





1 


\0 -i 


0) 




p~ 1 (z + b + r), x 





/n 

r 2 



Vo 



\ 






±W (z + b) 2 + p 2 , and r, is also of the form (14.9|) with a constant 



6^ 7^ 6. The conditions for U and V restrict these matrices to the form 



U 



ui u 3 
u 2 n 4 



-u 2 
Ul 



— M4 

U3 



V 



V3 «4 

-v 2 Vl 

\-V 4 V3J 



The latter condition may indeed be violated, as shown in Remark I5T2 
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with constants itj, «j. If r± = =■ r (i.e. b[ = b' 2 ), it turns out that g does not depend on r 
and Vi, and we obtain 

£ _ t + iaQ + fr)?- 1 - (m + it) 
~ t + ia(z + 6)r -1 + (m + il) ' 

with the parameters 

1=^7 — 1 — r, tn - 



2(liili4 — U2U3) ' 2(uiti4 — 7X2^3) ' U\U^ — U2U^ 

which satisfy m 2 — a 2 + ( 2 = 0. £ is the Ernst potential of an extreme Kerr-NUT space-time. 

Example 5.3 (Multi-Kerr-NUT). According to Remark 15 -H there is a simple way to superpose 
solutions by block-diagonally composing their matrix data. Let 



Pi = (v - 1 I i -Ri 



1 1>i 
-Uj 1 

r< 
-1/r, 

where pi ^ r^, i,k = 1, . . . , N , are given by (|4.9p . and either 6 j , 6 • , u j , n j E R or ^ = b. L G C, 
jj' = — jj, Uj = Uj G C (cf. Example 15. ip . Set 



U = (U 1 ,...,U N ), V 



\v N/ 



and P = block-diag(Pi, . . . , Pjv), R = block-diag(J?i, . . . , Rn), T = block-diag(7, . . . ,7). 
With go = I2, all assumptions of part (1) of Proposition 15.11 hold, hence with these data (I4.8P 
determines a family of solutions of the Ernst equation. Obviously, such a solution is a superpo- 
sition of N (non-extreme, respectively hyperextreme) Kerr-NUT solution^]. More generally, in 
the same way we can superpose any number of solutions with matrix data of the form given in 
Example 15.21 

5.2 Solutions of the Ernst equations in the Einstein— Maxwell case 

Choosing 

/0 i \ 

(5.15) 







i 





H 


i 















— Im-2 



and writing 




with a complex function £ and a complex (m — 2)-component vector <1>, (|5.2p takes the form 

1 §(£-£ + 2$t$) y/2& 



£+£+2&<$> 



13 See e.g. [6j [24] [51] for other derivations, and also [52], as well as the references cited there. 
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(also see [531 HZ])- We have 

1 - iff 2 i , T 1 - x 

where T denotes transposition. In the following we consider the case m = 3, where ()4.2p becomes 
the system of Ernst equations 

(Re£ + $<S>)(d 2 p + p~ l d p + 8 2 z )£ = (£ p ) 2 + (£ z ) 2 + 23>[<Z> p £ p + 

(Re 5 + $$) (d 2 + p- 1 ^ + «9 2 )$ = £ p $ p + + 2$ [($ p ) 2 + ($ z ) 2 ] , 

which determine solutions of the stationary axially symmetric Einstein-Maxwell equations (with- 
out further matter fields). If £ = 1 and $ = 0, then g reduces to 

A o o\ 

30= 1 , (5.16) 
\o -\j 

which corresponds to the Minkowski metric. 

Example 5.4 (Demiariski-Newman). Let n = 2 and 

-i o)' P= (o R= (o -1/r)' 

with p, r as in (I4.9P with constants 6, 6'. Solving gojU = UT and TV = V<?o7i an d recalling 
that U and V enter the solution formula (j4.8j) only up to an overall factor, leads to 

t/=|» 1 ], V-| » -;). (5.17) 

According to Proposition 15.11 part (1), in order to obtain solutions of the Ernst equations it 
remains to determine conditions under which g is Hermitian. By explicit evaluation one finds 
that this is so if one of the following sets of conditions is satisfied. 

(1) b,b' G R and 

U I 1 . Q . Q . Q > 

st = — 2— Imu, + Imu + |u — i lmv = 0, 2Imu+|s| =0. (5.18) 

u + i 

(2) V = b, j' = —j, v = u and t = s. 

Without restriction of generality we can set b' = —b, so that p and r are given by (|5.12|) . The 
Ernst potential £ is again of the form (|5.14|) . where now 

, 1 + uv — st 1 — UV U + V 

a = -b — — , l = jb — — , m = -jb- 




u — v + ist u — v + ist u — v + ist 

The second Ernst potential is given by 

$ = 2 ^ ± jgjg) (5 19) 

r + - jft- - if (r+ + jft-) + 2(m + it) ' 1 ' ' 



where 



jb s(v — i) + t(u + i) . jb s(v — i) — t(u + i) 



y/2 u — v + ist y/2 u — v + ist 
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In both cases, the parameters a, I, m, q e , q m are real and satisfy 

2 2 , r2 2 2 i2 



nr 



7 2 



Cases (1) and (2) correspond to a non-extreme, respectively hyperextreme, Demiariski-Newman 
space-time (see e.g. [50]). q e and g m are the electric and magnetic charge, respectively. Whe- 
reas (2) can be neatly expressed via (|5.7p . we have been unable so far to find a corresponding 
formulation for the conditions (1) in terms of the matrices P, R, U, V, also see Remark [57 



Example 5.5 (Harrison transformation). We can generate solutions of (|5.18p via a Harrison 
transformation [48 1 149]. A non-extreme Kerr-NUT solution (without charge) corresponds to the 
data 




Ur 



with real uq, vq. The matrix 

1 

H 




1 

-v 1 



iV2c 
V2c 



-1 - \c 

with c G C satisfies dSHD and (j5TT0]h Then 17' = 77t/ and V = satisfy g jU' = U'T 

and r V' = V'go-y, since ?7 o and Vq satisfy these conditions. Without effect on the solution of 
the chiral model, we can rescale these matrices to 



U 



V 



1 



1 + iuolcl 



-(n - i|c| 2 ) N 



1 + iuo|c| 2 
1 

1 — it>o|c| 2 



Uq — \\c\ 1 + mo c 



U )C 1 



>/2(i 
u + i|c 



1 — ivo|c|" 
-(v + i\c\ 2 ) l-if |c| 



2 



\/2(i + v )c 

-i\/2(i + vo )cy ' 



which have the form (j5.17j) and indeed satisfy (I5.18p . Using (I5.13p . the resulting Ernst poten- 
tials £ and $ can be written in the form (|5.14j) . respectively (|5,19p . where now 



1 + U V 

i . 

U - V 



l = jb 



1 - UqVq 1 + \c\ 
Uq-Vq 1 - Id 1 



m 



-jb 



U + V 1 + |c| 
Uq-Vq I - \c\ 



and 



2j^[(nof o — 1) Re c — (uq + f o) Im c] 
(1 - |c| 2 )(n - vq) 
2jb[(uQ + fo) Re c + (uqVq — 1) Im c] 

(1 - \c\ 2 )(uq - Vq) 

Example 5.6 (hyperextreme multi-Demiahski-Newman) . Let 

\ P = block-diag(Pi, ...,P N ) 

R = block-diag(J?i, . . . , Rn) 
U = {U U ...,U N ), 



/0 i 
-i 



\ 



i 
-i 0/ 



\v N ) 
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with 




-Ui 







) 



) 



Ri — 



( 



( 



-Ui 







1 



1 isi 



i/n 



o 



) 



) 



where r% = —1/pi- Then the conditions of Proposition 15. II are obviously satisfied (with 7 and go 
given by (|5.15p . respectively (|5.16p ). It follows that fj4.8|) determines a solution of the (Einstein- 
Maxwell-) Ernst equations. This is a superposition of N hyperextreme Demiahski-Newman 
solutions. 

Remark 5.4. Whereas the hyperextreme mti/ii-Demiahski-Newman solutions are obtained in 
a straightforward way, this is not so in the non-extreme case. So far a suitable condition on the 
matrix data is lacking. Similar problems are known in other approaches, see e.g. [SJ Section 34.8]. 

Example 5.7. Let n = 2N. With the choices made in Remark I5.2[ Proposition 15.11 part (1), 
implies that the expression given by (|4.8p with g$ in (|5.16p satisfies (jg) 2 = I3 and tr(7<jf) = 1. 
In order to obtain a solution of the Ernst equations, it suffices to arrange that g is Hermitian. 
A sufficient condition is given by part (2) of Proposition 15.11 This leads to a huge family of 
solutions of the Einstein-Maxwell equations. The hyperextreme Demiahski-Newman solution 
is just the simplest example in this family. Furthermore, such solutions can be superposed in 
the simple way described in Remark 15. II and Example 15.61 An exploration of the corresponding 
space-times would be a difficult task. 

6 Conclusions 

We addressed the m x m non- autonomous chiral model in a new way, starting from a very 
simple and universal solution generating result within the bidifferential calculus approach. This 
resulted in an infinite family of exact solutions for any matrix size m, parametrized by matrices 
subject to a Sylvester equation. To solve the latter is a well-studied and fairly simple problem. 
At least in the compact form presented in this work, according to our knowledge these solutions 
have not appeared previously in the literature. 

The non-autonomous chiral model originally appeared in reductions of Einstein's equations. 
We demonstrated in Section [5] that the "multi-solitons" on a flat background, known in the case 
of stationarity, axial symmetry and vacuum, respectively electrovacuum, are indeed contained 
in the family of solutions that we obtained in Section [J] for the non-autonomous chiral model 
equation. More precisely, we found conditions to be imposed on the matrices that parametrize 
the latter solutions such that (in the cases m = 2, respectively m = 3) they become solutions of 
the Ernst equation(s) of general relativity. Only in the case of non-extreme multi-Demiariski- 
Newman solutions we were not (yet) able to find a corresponding characterization of the matrix 
data. 

Beyond the solutions found e.g. by Belinski and Zakharov, which in the present work corre- 
spond to diagonal matrices P and R, there are solutions associated with non-diagonal matrix 
data. It may well be that such solutions can be obtained alternatively e.g. in the Belinski- 
Zakharov framework with a dressing matrix involving higher order poles, or by taking suitable 
limits where some poles coincide. In any case, our approach yields these solutions directly. 
Moreover, relaxing the spectral condition for the matrices P and R, the Sylvester equation has 
further solutions, provided that the matrix VU on its right hand side is appropriately chosen. 
This is another possibility to obtain new solutions. Finally, we should mention the possibility to 
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make sense of the limho n — > oo (where nxnis the size of the matrices that parametrize the 
solutions). In conclusion, at present it is not quite clear what the generated class of solutions 
really embraces. 

Moreover, using the original method of Belinski and Zakharov, in the Einstein-Maxwell case 
no appropriate reduction conditions could be found (cf. [SI]), and a different approach had to 
be developed [MIE]- We had less problems in this respect. 

On the other hand, the Belinski-Zakharov approach, the modified approach of Alekseev [23] 
in the Einstein-Maxwell case, and others can also be used to generate "solitons" on a non- 
flat background. Perhaps a corresponding extension of Proposition 13.11 exists. This is also 
suggested by the relation with Darboux transformations in Appendix [A] In any case, here we 
have a limitation of Proposition 13.11 (which is not a limitation of the bidifferential calculus 
framework, which offers various methods pQ), but we have the advantage of a very simple and 
general result that covers physically interesting cases. 

The appearance of a Sylvester equation is a generic feature of the solution generating result 
formulated in Proposition 13.11 and in Appendix IA1 (also see [1]). Sylvester equations and their 
simplest solutions, Cauchy-like matrices, frequently appeared in the integrable systems liter- 
ature. But this is the first time we came across a Sylvester equation involving non-constant 
matrix data. A particularly nice feature is the fact that solutions can be superposed by simply 
composing their matrix data into bigger block-diagonal matrices. The corresponding Sylvester 
equation still has to be solved, but a unique solution exists if we impose a not very restrictive 
spectral condition on these matrix data. 

In Appendix[B]we recovered two familiar Lax pairs for the non-autonomous chiral model from 
the general linear equation (|2.7p in the bidifferential calculus framework. Our way toward exact 
solutions in Section [4] is more closely related to Maison's Lax pair than to that of Belinski and 
Zakharov. We eliminated the ^-dependence, whereas in the Lax pair of Belinski and Zakharov 
the 6>-dependence is kept and it involves derivatives with respect to this "spectral parameter" . 

Our results extend beyond the Einstein-Maxwell case and are also applicable to higher- 
dimensional gravity theories (see e.g. HJH HH H2J H31 E31 H5J EEl EBl US] ) ■ Besides that, other 
reductions of the non-autonomous chiral model (for some m) are of interest, see e.g. [201121]. and 
the set of solutions that we obtained in this work will typically be reducible to solutions of them. 

Since Proposition 13.11 actually generates solutions of a Miura transformation equation, we 
obtained simultaneously solutions of the Miura-dual of the non-autonomous chiral model equa- 
tion: under the assumptions of Proposition 14.11 <j) = UX V solves the equation (|4.3p . The 
further results in Section [H in particular Example 14.11 then provide us with explicit families of 
solutions. 

Proposition 13. 1[ respectively Theorem I A. 11 can actually be formulated and proved without 
explicit use of the two nonlinear equations involving only P, respectively R. In such a for- 
mulation, the theorem generates solutions of the nonlinear integrable equation (12.3p . respec- 
tively (|2.4p . from solutions of linear equations. However, the equations for P and R arise as 
integrability conditions of the latter. In previous work [TJ [32J [33J [34] , we chose P and R as d- and 
d- constant matrices, which indeed reduces the equations that have to be solved to only linear 
ones, and we recovered (and somewhat generalized) known soliton solution families. In case of 
the non-autonomous chiral model and, more specifically, its reduction to the Ernst equation, it 
turned out to be necessary to go beyond this level, and thus to consider genuine solutions of 
the nonlinear equations for P and R, in order to obtain relevant solutions like those associated 
with multi-Kerr-NUT space-times and their (electrically and magnetically) charged generaliza- 
tions. This also suggests a corresponding application of the theorem (or Proposition 13. It) to 
other integrable PDDEs. 



See e.g. [55] [56] for results on the operator Sylvester equation. 
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A Via a Darboux transformation and a projection 
to a non-iterative solution generating result 

Lemma A.l. Let P be invertible. The transformation 

(4>,g) ^ g')=(<t> + XPX~\XPX- 1 g), (A.l) 

where X is an invertible solution of (j2.7|) with A = dip = (dg)g~ 1 , and dP = (dP)P, maps 
a solution of the Miura transformation equation (|2.10|) into another solution. 

Proof. Using (|2.7[) and dP = (dP)P, a direct computation leads to 

(dg^g'" 1 -d(/)' = A-d(p-XPX~ 1 [A-(dg)g~ 1 ]XP~ 1 X~ 1 , 

which vanishes if A = d(f> = (dg)g~ 1 . ■ 

(jA.ip is an essential part of a Darboux transformation, cf. [Tj. In the following we will use 
this result to derive a theorem which essentially reduces to Proposition 13.11 see Remark IA.21 

Lemma A. 2. Let ((f), g) be a solution of the Miura transformation equation (|2.10p in 
M&t(n,n,B). Let U € Mat(m, n,B) and V G Mat(n,m, B) be d- and d- constant. If 

<f> = VU<f>, (A.2) 
with some <f> 6 Mat(n, n, B), then 

4> = U$V, g =(Ug- 1 Vy 1 (A.3) 
solve the Miura transformation equation (|2,10p in Mat(m,m,£>). 
Proof. Since (4>,g) is assumed to solve (|2.10p . we have 

dfT 1 = -g- l d<t> = -g- l VUd$. 
Multiplying by U from the left and by V from the right, we obtain 

dcT 1 = -g^dcp, 

which is equivalent to (|2.10|) . ■ 

Theorem A.l. Let (—R,S) be a solution of the Miura transformation equation (|2.10p in 
Mat(n,ra,£>), i.e. 

dS = -{dR)S, (A.4) 

and S invertible. Let X be an invertible solution of the linear equation (|2.7p . now in Mat(n, n, B) 
and with invertible P, hence 

dX = (dX)P - (dR)X, dP = (dP)P. (A.5) 

In addition we require that 

XP - RX = VUY, dR = RdR, dY = (dY)P, (A.6) 

where U E Mat(m,n, B) and V G Mat(n, m,B) are d- and d-constant, and Y G Mat(n, n,B). 
Then also 

cb = UYX~ 1 V and g = (US^XP^X^Vy 1 (A.7) 
solve the Miura transformation equation (|2.10p . and thus (|2,3p . respectively ([22 
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Proof. Since we assume that (— R, S) solves the Miura transformation equation (|2.10p in 
Mat(n,n,B), according to Lemma lA.ll this also holds for the pair 

<f} = -R + XPX~ 1 , g = XPX 1 S. 

Using the first of we find that (|A~2j) holds with cj> = YX~ l . Now (jAT3j) yields the 

asserted formulas for (j) and g. According to Lemma [A. 2\ <p and g solve the Miura transformation 
equation (|2.10p . 

Together with (|A.5|) . the first of (|A.6P implies 

VU(dY - (dY)P) = (RdR-dR)X, 

which is satisfied if the last two conditions of (|A.6P hold. ■ 

Remark A.l. This theorem generalizes a previous result in [I], which has been applied in 
[H [32"1 [33l 133] with d- and d-constant P, R, in which case only linear equations have to be 
solved in order to generate solutions of (|2.3p . respectively (|2.4p . 

The above derivation shows that the theorem may be regarded as a combination of a Dar- 
boux transformation (Lemma lA.ip . on the level of matrices of arbitrary size, and a projection 
mechanism (Lemma IA.2p . The projection idea can be traced back to work of Marchenko |57j . 
More generally, the above result can be formulated in terms of suitable operators, replacing the 
matrices that involve a size n. 

The next remark shows that, with mild additional assumptions, Theorem I A. II reduces to 
Proposition 13.11 

Remark A. 2. The transformation 

X t-> XQ. Y^YQ, P^Q^PQ, 

with any invertible Q £ M.&t(n,n,B), leaves the expressions for (j) and g in (|A.7p invariant. 
This is then also a symmetry transformation of (|A.5P and (|A.6P if dQ _1 = (dQ~ 1 )P. As 
a consequence, under the assumptions that Y is invertible, without restriction of generality we 
can set Y = I, where / is the n x n identity matrix. Then (j) is given by the expression in (|3.2p . 

We further note that (|A.4p and the second of (|A.6P imply d(RS) = 0. Assuming that R 
is invertible, we thus have S = i? _1 C with an invertible d-constant C. The expression for g 
in (|A.7p now takes the form 

g = (uc^Rxixpy^y 1 = (i - u(xpy 1 vy l (uc^vy 1 , 

assuming temporarily inver tibility of UC^V. Together with cf), this remains a solution of (]2. lQf) 
if we drop the last factor, so that 

g = {i-u(xpy 1 vy 1 . 

This expression also makes sense without the above additional invertibility assumptions. We 
can still translate it into a simpler form. From the first of (|A.6[) . which now has the form of the 
last of (|3.ip . we obtain 

(Rxy 1 - (xpy 1 = (Rxy 1 vu(xpy 1 . 

Multiplication by U from the left and by V from the right, and use in our last formula for g, 
leads to the expression for g in (|3.2p . 
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B Linear systems for the non-autonomous chiral model 

B.l Maison's Lax pair 

Using the bidifferential calculus determined by (|4.1|) . (|2.7|) with A = (dg)g~ 1 takes the form 
X p + p- x X B = {g p + p- 1 g e )g- 1 X - X z Pe 9 , X z = g z g' l X + (X p - p~ 1 X )Pe e , 

and (|2.8|) reads 

e~ e (P p + p~ l P e ) = -P Z P, P z = e e (P p - p~ 1 P e )P 

Disregarding a constant solution (cf. Section |B.2|) . we can eliminate the 0-dependence in the 
latter equations via 

P = e- e P, 
with P independent of 9, and obtain 

P p = p~ l P{l - P 2 ) (I + P 2 )~\ P z = 2p~ x P 2 {l + P 2 ) -1 . 
Furthermore, setting 

X = e^ e X, g = e c > e g, 

with X , g independent of 9, the above linear system becomes 

X p (l + P 2 ) = (~g p + c 2 p- 1 ~g)~g- 1 X-~g z ~g- 1 XP-c 1 p- 1 X(l-P 2 ), 
X Z (I + P 2 ) = g.g^X + (g p + c 2 p- l ~g)~g~ l XP - 2c 1 p" 1 XP. 

Choosing 

P = pl 

with a function p(p, z), the equations for P can easily be integrated, which results irj^l 

p = P" 1 (z + b ± V(2 + fe) 2 +p 2 ) , 
where b is an arbitrary constant. In terms of 

X = ~g- l X 

the above linear system, simplified by setting c\ = c 2 = 0, then takes the form 

X p = - P n (ff -1 fe +pg~ 1 g p )X, X z = - P (g- l g p -pg- l g z )X. 
1 + p z 1 + p z 

This system is equivalent to a linear system for the non-autonomous chiral model, first found 
by Maison in 1979 0] (also see [21]). 



3 We note that p is Pi in Section 1472 
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B.2 The Belinski-Zakharov Lax pair 

Using instead of 6 the variable 

A = -pe e , 
(|4.ip translates into 

d/ = -Mi - P' 1 XM2, d/ = -{pX^fp + 2/ A )Ci + Mi- 
We consider the linear system (|2.7|) with P = I, which trivially solves (|2.8|) . i.e. 

dX = AX + dX 
Writing 

the integrability condition ()2.9p takes the form 

Bp - A, = [A, B], (pA) p + ( P B) Z = 0, 
assuming that A, B are A-independent. Solving the first (zero curvature) condition by 

A = g p g~ 1 , B = g z g~ 1 , 
the second equation becomes the non-autonomous chiral model equation 

{pg P g~ l ) p + {P9z9~ l ) z = 0. 
The above linear equation leads to 



pU + W 2p\ P V - XU 2A 2 

p 2 + A 2 p 2 + A 2 A ' z p 2 + X 2 + p 2 + X 2 ' 



where 

W = P A, V = pB. 

This is the Belinski-Zakharov Lax pair [6] (also see Chapter 8]). We note that the "spectral 
parameter" A has its origin in a coordinate of the self-dual Yang-Mills equation. We also note 
that A = (dg)g~ l (using g x = 0). 

C Some proofs 

C.l Proof of Lemma 14.11 

~ 2 

(1) Assuming that I + P is invertible, the system (|4.4p can be decoupled into 

p p = P - l P{i - P 2 ) {i + py\ P z = 2 P - l P\i + P 2 y\ (c.i) 

which can also be written as 

(p-\ = - p - x p-\i - p 2 ) (i + b 2 )- 1 , (p~ l ) z = -ip~ i (i + p 2 r\ 
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assuming that P is invertible. Subtraction yields 

(P - P~\ = -p" 1 (P - P" 1 ), (P - P~\ = 2p~% 
which can be integrated to 

P- P" 1 = 2p~ l (zI + B), (C.2) 

with a constant matrix B. This implies (|4.5|) . 

(2) Let P satisfy (|4.5p with a constant matrix P. Then P is invertible, since the existence 
of a non- vanishing vector annihilated by P would be in conflict with (|4,5p . Thus (|C.2p holds, 
which implies [P, P] = 0. Differentiation of (|4.5p with respect to p, and elimination of z J + B 
with the help of (|4.5p or equivalently (|C2[) . leads to 

= P P P + PP P + 2/)- 2 (z/ + B)P - 2p~ 1 (zl + P)P P = (P + P _1 )P P + P~H p2 " J ) . 

~ ~ 2 

where we used the assumption [P p , P] = 0. If J + P is invertible, the resulting equation is the 
first of (jC.ip . In the same way we obtain the second of (jC.ip . (jC.ip is equivalent to (|4.4p . 

C.2 Proof of Proposition 15.11 

Using ([477]) . ([53D and ([53]) . we find 

r(xp + rpxr) - (xP + rPxr)P = o, 

so that the spectrum condition implies 

TRXT = -XP. (C.3) 
With the help of this result we obtain 

m (l + UiRX^V) = g Ql + UT(RXy l V = g Ql - UiTXP^V 

= (i-uixPy^gor 

Using (go"/) 2 = I, the condition (jg) 2 = I for (]4.8f) is therefore equivalent to 

(/ - utxPyW) (i + [/(Mj^y) = /. 

Expanding the left hand side and using the Sylvester equation (|4.7|) to eliminate VU, this indeed 
turns out to be satisfied. To complete the proof of (1), it remains to derive the trace formula. 
Using ([P]) . (ET7D and (i03l) . we obtain 

tr( 7 $) - tr( 7 5o) = ^((M)- 1 Vg 0l U) = ^((M)- 1 Vl/r) 

= tr((PX)~ 1 (XP - M)r) = -tr(r) + tr((Px)~ 1 xPr) 
= -tr(r) - tT^Rx^rhx) = -2tr(r). 

In order to prove (2), we consider the Hermitian conjugate of the Sylvester equation (|4.7p . 
By use of (I5.7p . and with the help of gojU = UT, TV = Vgo"/ and (go"/) 2 = I, it takes the 
form 

vu = -(rx t r)P + p(rx t r). 

By comparison with the original Sylvester equation, the spectrum condition allows us to conclude 
that 

X ] = -TXT. 

Together with (1C.3P this implies 

(pjf ) f = x j k^ = -rxr 2 pr = -rxPr = kx. 

It follows that U(RX)~ 1 Vgo is Hermitian, and thus also g given by (|4.8[) . 



The Non-Autonomous Chiral Model in the Bidifferential Calculus Framework 



25 



Acknowledgements 

We would like to thank Vladimir S. Manko and anonymous referees for helpful comments. 
During the course of this work, N.K. has been at the Max-Planck-Institute for Dynamics and 
Self-Organization in Gottingen. 

References 

Dimakis A., Miiller-Hoissen F., Bidifferential graded algebras and integrable systems, 
\Discrete Contin. Dyn. ^1) 2009 (2009), suppl., 208-219, lajXiv:0805.4553l 

Kinnersley W., Symmetries of the stationary Einstein-Maxwell field equations. I, \J. Ma th. Phys. 18 (1977), 
1529-1537. 

Maison D., Are the stationary, axially symmetric Einstein equations completely integrable?, Phys. Rev. Lett. 
41 (1978), 521-522. 

Maison D., On the complete integrability of the stationary, axially symmetric Einstein equations, 
| J. Math. PhyT\ 20 (1979), 871-877. 

Belinskil V., Zakharov V., Integration of the Einstein equations by means of the inverse scattering problem 
technique and construction of exact soliton solutions, Soviet Phys. JETP 48 (1978), 985-994. 

Belinskil V., Sakharov V., Stationary gravitational solitons with axial symmetry, Soviet Phys. JETP 50 
(1979), 1-9. 

Belinski V., Verdaguer E., Gravitational solitons, Cambridge Monographs on Mathematical Physics, Cam- 
bridge University Press, Cambridge, 2001. 

Stephani H., Kramer D., MacCallum M., Hoenselaers C, Hcrlt E., Exact solutions of Einstein's field equa- 
tions, Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge, 2003. 

Lee S.-C, Generalized Neugebauer-Kramer transformations for nonlinear sigma models, \Phys. Lett. B\ 164 
(1985), 75-79. 

Lee S.-C, Soliton transformations for axially symmetric higher-dimensional gravity. II. Belinkii-Zakharov 
A-soliton transformations, | J. Math. Phys.\ 28 (1987), 901-904. 

Wei Y.-H., Extended Ernst equation and five-dimensional axisymmetric Kaluza-Klein solutions, 
\Classical Quantum Gravity] l8 (2001), 2163-2170. 

Harmark T., Stationary and axisymmetric solutions of higher-dimensional general relativity, \Phys. Rev. D\ 
70 (2004), 124002, 25 pages, [hepTh/0408141| 



Pomeransky A. A., Complete integrability of higher-dimensional Einstein equations with additional symme- 
try and rotating black holes, \Phys. Rev. J>| 73 (2006), 044004, 5 pages, |hep-th/0507250| 

Azuma T., Koikawa T., An infinite number of stationary soliton solutions to the five-dimensional vacuum 



Einstein equation, \Progr. Theoret. Phys.\ 116 (2006), 319-328, [hep-th/0512350 



Iguchi H., Mishima T., Solitonic generation of vacuum solutions in five-dimensional general relativity, 



\Phys. Rev. D\ 74 (2006), 024029, 17 pages, [hep^t h/0605090 



Tomizawa S., Morisawa Y., Yasui Y., Vacuum solutions of five dimensional Einstein equations generated by 



inverse scattering method, \Fhys. Rev. D] 73 (2006), 064009, 8 pages, hep-th/0512252 

Yazadjiev S.S., 5D Einstein-Maxwell solitons and concentric rotating dipole black rings, Phys. Rev. D 78 
(2008), 064032, 11 pages, larXiv:0805.1600l 

Emparan R., Reall H.S., Black holes in higher dimensions, Living Rev. Relativity 11 (2008), 6, 87 pages, 
larXiv:0801. 34711 

Figueras P., Jamsin E., Rocha J.V., Virmani A., Integrability of five-dimensional minimal supergravity and 
charged rotating black holes, Classical Quantum Gravity 27 (2010), 135011, 37 pages, arXiv:0912.3199 

Mikhailov A.V., Yaremchuk A.I., Cylindrically symmetric solutions of the non-linear chiral field model (<r 
model), \Nuclear Phys. B 202 (1982), 508-522. 

Gutshabash E.Sh., Lipovskii V.D., Exact solutions of a nonlinear sigma-model in a curved space and the 
theory of magnetically ordered media with variable saturation magnetization, J. Math. Sci. 77 (1995), 
3063-3068. 

Alekseev G.A., Af-soliton solutions of Einstein-Maxwell equations, JETP Lett. 32 (1980), 277-279. 



26 



A. Dimakis, N. Kanning and F. Miiller-Hoissen 



[23] Alekseev G.A., Exact solutions in the general theory of relativity, Proc. Steklov Inst. Math. 176 (1988), 
215-262. 

[24] Neugebauer G., A general integral of the axially symmetric Einstein equations, J. Phys. A: Ma th. Gen.\ lS 
(1980), L19-L21. 

[25] Neugebauer G., Kramer D., Einstein-Maxwell solitons, \L Phys. A: Math. Gen.\ 16 (1983), 1927-1936. 

[26] Kramer D., Neugebauer G., Matos T., Backlund transforms of chiral fields, \J. Math. Phys.\ 32 (1991), 
2727-2730. 

[27] Korotkin D., Finite-gap solutions of the stationary axisymmetric Einstein equation in vacuum, 
| Theoret. and Math. Phys] 77 (1989), 1018-1031. 

[28] Manko V., Martin J., Ruiz E., Extended family of the electrovac two-solitons for the Einstein-Maxwell 
equations, \Phys. Rev. J>] 51 (1995), 4187-4191. 

[29] Ruiz E., Manko V., Martin J., Extended iV-soliton solution of the Einstein-Maxwell equations, \Phys. Rev. D\ 
51 (1995), 4192-4197. 

[30] Masuda T., Sasa N., Fukuyama T., Neugebauer-Kramer solutions of the Ernst equation in Hirota's direct 
method, ! J. Ph V s - A: Math. Gen.\ 31 (1998), 5717-5731. 

[31] Burtsev S.P., Zakharov V.E., Mikhailov A.V., Inverse scattering method with variable spectral parameter, 
| Theoret. and Math. Phys] 70 (1987), 227-240. 

[32] Dimakis A., Miiller-Hoissen F., Solutions of matrix NLS systems and their discretizations: a unified treat- 
ment. ITrwerse Problem^ 26 (2010), 095007, 55 pages. 

[33] Dimakis A., Miiller-Hoissen F., Bidifferential calculus approach to AKNS hierarchies and their solutions, 
\SIGMA\ 6 (2010), 055, 27 pages, larXrv: 1004. 16271 

[34] Dimakis A., Kanning N., Miiller-Hoissen F., Bidifferential calculus, matrix SIT and sine-Gordon equations, 
Acta Polytechnica 51 (2011), 33-37. [arXivilOll. 17371 

[35] Witten L., Static axially symmetric solutions of self-dual SU(2) gauge fields in Euclidean four-dimensional 
space, \Ph~ys. Rev. J>| 19 (1979), 718-720. 

[36] Forgacs P., Horvath Z., Palla L., Generating the Bogomolny-Prasad-Sommerfield one-monopole solution 
by a Backlund transformation, \Phys. Rev. Lett . 45 (1980), 505-508. 

[37] Bais F.A., Sasaki R., On the complete integrability of the static axially symmetric self-dual gauge field 
equations for an arbitrary group, Nuclear Phys. B 195 (1982), 522-540. 

[38] Ward R.S., Stationary axisymmetric space-times: a new approach, Gen. Relativity Gravitation 15 (1983), 
105-109. 

[39] Mason L.J., Woodhouse N.M.J., Integrability, self-duality, and twistor theory, London Mathematical Society 
Monographs. New Series, Vol. 15, Oxford Science Publications, The Clarendon Press, Oxford University 
Press, New York, 1996. 

[40] Klein C, Richter O., Ernst equation and Riemann surfaces. Analytical and numerical methods, Lecture 
Notes in Physics, Vol. 685, Springer- Verlag, Berlin, 2005. 

[41] Kanning N., Integrable Systeme in der Allgemeinen Relativitatstheorie: ein Bidifferentialkalkiil-Zugang, 
Diploma Thesis, University of Gottingen, 2010. 

[42] Higham N.J., Computing real square roots of a real matrix, Linear Algebra Appl. 88/89 (1987), 405-430. 

[43] Bernstein D.S., Matrix mathematics. Theory, facts, and formulas, 2nd ed., Princeton University Press, 
Princeton, NJ, 2009. 

[44] Zakharov V.E., Mikhailov A.V., Relativistically invariant two-dimensional models of field theory which are 
integrable by means of the inverse scattering problem method, Soviet Phys. JETP 47 (1978), 1017-1027. 

[45] Eichenherr H., Forger M., More about non-linear sigma models on symmetric spaces, Nuclear P hys. j?| 164 
(1980), 528-535. 

[46] Mazur P.O., Proof of uniqueness of the Kerr-Newman black hole solution, \J. Phys. A: Math. Gen. 15 (1982), 
3173-3180. 

[47] Giirses M., Inverse scattering, differential geometry, Einstein-Maxwell solitons and one soliton Backlund 
transformations, in Solutions of Einstein's Equations: Techniques and Results (Retzbach, 1983), Editors 
C. Hoenselaers and W. Dietz, \Lecture Notes in P hys., Vol. 205, Springer, Berlin, 1984, 199-234. 

[48] Harrison B.K., New solutions of the Einstein-Maxwell equations from old, | J. Math. Phys] 9 (1968), 1744- 
1752. 



The Non-Autonomous Chiral Model in the Bidifferential Calculus Framework 



27 



Kinnersley W., Generation of stationary Einstein-Maxwell fields, [J. Math. P hys. 14 (1973), 651-653. 

Aguilar-Sanchez J., Garcia A. A., Manko V.S., Demiariski-Newman solution revisited, Gravit. Cosmol. 7 
(2001), 149-152, [gr7qc/0106011| 



Kramer D., Neugebauer G., The superposition of two Kerr solutions, Phys . Lett. A\ 75 (1980), 259-261. 

Neugebauer G., Hennig J., Non-existence of stationary two-black-hole configurations, 
| Gen. Relativity Gravitatib~h] 41 (2009), 2113-2130. rarXiv:0905.4179l 

Giirses M., Xanthopoulos B., Axially symmetric, static self-dual SU(3) gauge fields and stationary Einstein- 
Maxwell metrics, \Phys. Rev. D\ 26 (1982), 1912-1915. 

Alekseev G., Thirty years of studies of integrable reductions of Einstein's field equations, arXiv:1011.3846 

Bhatia R., Rosenthal P., How and why to solve the operator equation AX — XB — Y, 
\Bull. London Math. 5oc~l 29 (1997), 1-21. 

Garimella R.V., Hrynkiv V., Sourour A.R., A solution of an operator equation related to the KdV equation, 
\Linear Algebra Appl.\ 418 (2006), 788-792. 

Marchenko V., Nonlinear equations and operator algebras, Mathematics and its Applications, Vol. 17, D. Rei- 
del Publishing Co., Dordrecht, 1988. 



